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ABSTRACT 


The seismic microzonation of the urban area of Ambato, was made in 2018, 
in a probabilistic and a deterministic way. This last form of calculation is the one 
presented in the first part of the article; for this purpose, 3 geological faults and 3 
strong motion equations were considered. 


Therefore, thinking of incorporating the rural area than is approximately four 
times than the urban area, 4 more geological faults and 3 new attenuation models 
have been incorporated. In such a way that the Methodology of calculation of the 
seismic hazard in a deterministic way of the urban and rural area of Ambato 
considers 7 geological faults and 6 strong motion equations. It is a preliminary work 
that has the character of an exploratory since it requires studies of geophysics, 
geotechnics, and geology in rural areas. 

Keywords: Seismic hazard in a deterministic way, Strong motion equations, Control 
spectra, Horizontal component 


METODOLOGIA PARA EVALUAR LA PELIGROSIDAD 
SÍSMICA DEL ÁREA URBANA Y RURAL DE LA 
CIUDAD DE AMBATO EN FORMA DETERMINÍSTICA 


RESUMEN 


La microzonificación sísmica del área urbana de la ciudad de Ambato fue 
realizada en 2018, en forma probabilística y en forma determinística. Esta última 
forma de cálculo es la que se presenta en la primera parte del artículo; para el efecto 
se consideraron 3 fallas geológicas y 3 ecuaciones de movimientos fuertes. 


Por lo tanto , pensando en incorporar el área rural que aproximadamente es 
cuatro veces el área urbana, se han incorporado 4 fallas geológicas más y 3 nuevos 
modelos de atenuación. De tal manera que la Metodología de cálculo de la 
peligrosidad sísmica en forma determinística del área urbana y rural de Ambato 
considera 7 fallas geológicas y 6 ecuaciones de movimientos fuertes. Es un trabajo 
preliminar que tiene el carácter de exploratorio ya que hacen falta estudios de 
geofísica, geotécnica y geológicos en el área rural. 
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Palabras Claves: Peligrosidad sismica en forma deterministica. Ecuaciones de 
movimientos fuertes, espectros de control, componente horizontal 


1. INTRODUCCION 


La mejor forma de mitigacion de los desastres naturales es proporcionar al 
Proyectista Estructural informacion confiable sobre los espectros que debe 
considerar en el diseno sismico de las estructuras, la misma que sea obtenida en 
base a la mejor informacion: geologica y sismica existente a la fecha (2020). Un 
factor fundamental es que en los modelos de atenuacion de los movimientos del 
suelo también conocidos como ecuaciones de movimientos fuertes estén más 
apegados al entorno en el cual se halla la ciudad de Ambato, como el ambiente 
tectónico sea similar al sitio de estudio o el rango de aplicabilidad el modelo sea lo 
suficientemente amplio. 


En este contexto, se debe manifestar que Ambato, se encuentra cerca de 
los volcanes Tungurahua y Cotopaxi. Por lo tanto, se debe utilizar un modelo de 
atenuación de movimiento del suelo, que tome en cuenta esta variable. La mayor 
parte de modelos existentes, no toman en cuenta la presencia de los volcanes, en 
el cálculo de la ley de atenuación o de los espectros de respuesta. 


Un modelo que si considera la presencia de los volcanes en el estudio de 
peligrosidad sísmica es el de Zhao et al. (2016) que será considerado en el estudio 
para determinar espectros ante la componente horizontal de movimiento del suelo. 


Pero no se debe pensar que los estudios de peligrosidad sísmica van a 
servir únicamente a las nuevas construcciones, nada más falso, los estudios sirven 
fundamentalmente para evaluar la vulnerabilidad sísmica de las estructuras 
existentes. 


Leyenda 
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Figura 1 Estudios de Peligrosidad sismica del area urbana de la ciudad de Ambato; 
a) Períodos de vibración del suelo; b) Velocidad de la onda de corte en los 
30 primeros metros. Fuente: Bonifaz y Haro (2018). 


Con esta introducción a continuación se indica el estado del arte sobre 
peligrosidad sísmica, en el área urbana de la ciudad de Ambato y para ello se debe 
indicar los trabajos realizados por Aguiar y Rivas (2018), destacando que estos 
fueron respaldados por una serie de estudios que se indican a continuación y cuyo 
detalle de cálculo están inmersos en la referencia indicada. 


En la figura 1 a, se muestra las curvas de igual período de vibración del 
suelo, destacando que estos varían desde 0.1 s, hasta 0.45 s., por lo que, desde el 
punto de vista estructural, son buenos suelos. Para que se entienda mejor se debe 
indicar que en la ciudad de Guayaquil los períodos de vibración del suelo llegan a 4 
segundos, donde la amplificación de las ondas sísmicas será considerable. 


En la figura 1 b, se muestran las curvas con igual velocidad de la onda de 
corte en los primeros 30 metros V;3/, que varían entre 300 y 450 m/s, que se pueden 
catalogar como buenos suelos. Se debe indicar que existen puntos muy localizados 
donde V,39 es muy bajo que corresponden a rellenos de quebradas, pero no son 
representativos, con respecto al área de la ciudad de Ambato, por lo que no se ha 
colocado un valor menor de V;z0. 


Para poder clasificar los perfiles de suelo sísmico de acuerdo con la Norma 
Ecuatoriana de la Construcción NEC-15 se necesita conocer otras variables, la 
principal el número de golpes del ensayo de penetración estándar Nso, los mismos 
que fueron realizados por AOC Ingeniería (2018); se determinó también el esfuerzo 
de corte del suelo, el ángulo de fricción interna, la cohesión de los suelos y las 
propiedades mecánicas del suelo. En base a toda esa información y a la indicada 
en la figura 1, se obtuvo la clasificación de los suelos que se presenta en la figura 
2. 

Con los periodos de vibración, encontrados en la ciudad de Ambato se ha 
elaborado curvas que tienen igual periodo, presentadas en la figura 1a. La mayor 
parte de la ciudad tiene periodos del suelo inferior a 0.27 segundos, además se 


presenta que todos los suelos tienen velocidades V;3ọ9, mayores a 300 m/s, lo que 
significa desde el punto de vista de la ingeniería sísmica que se trata de suelos 
buenos. 


Con el fin de analizar la zona rural de la provincia se clasifico a los suelos 
tipo “C”, y suelo tipo “D” para determinar los espectros de control, mediante un 
estudio de peligrosidad sísmica Determinístico; para el efecto en cada perfil de 
suelo, se construyó una malla de puntos equidistantes cada 500 metros y se clasificó 
de acuerdo al tipo de suelo, se halló espectros en cada punto, se considera que una 
de las tres fallas genera un sismo de magnitud máxima, se analizan sus valores y 
se analiza que falla genera las mayores ordenadas espectrales, determinando 
espectros promedios en cada perfil de suelo. 
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Figura 2 Clasificación de los suelos de la ciudad de Ambato desde el punto de vista 
de la Peligrosidad Sísmica. Fuente: Aguiar y Rivas (2018). 


2. SELECCIÓN DE FALLAS ACTIVAS 


Existen dos metodologías para evaluar la peligrosidad sísmica de una 
región, la una denominada probabilística en que se trabajan con zonas fuentes y 
la otra llamada deterministica en que se trabaja con cierto número de fallas 
geológicas. En este trabajo, el cálculo de la peligrosidad sísmica se lo hace en forma 
determinística utilizando para el efecto las fallas corticales. 


Se decidió trabajar con las fallas corticales, a la vista de los resultados 
encontrados en la desagregación de la peligrosidad sísmica para la ciudad de 
Ambato que se presenta en la figura 3, para el PGA y para un período de 1 s., para 
un período de retorno PR de 475 años, en la cual se observa que para el PGA son 
importantes las fallas corticales y son más importantes aún para otros períodos de 


vibración. 





PGA - PR 475 años SA (1.0 s) - PR 475 años 
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S € S 156-02 € 
O 106-02 2 O 2 
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Figura 3 Resultados de la desagregacion de la Peligrosidad Sismica para la ciudad 
de Ambato, para un periodo de retorno PR de 475 anos; a) Valores para 
el PGA; b) Valores para un periodo de 1 s. Fuente: Aguiar y Rivas (2018). 


Al observar la figura 3 a, se podría pensar en considerar también los sismos 
de subducción (Interfase) ya que la contribución inclusive es mayor que la 
proveniente de sismos corticales, pero para períodos diferentes de cero como el 
mostrado en la figura 3 b, la contribución notable al cálculo de la peligrosidad 
sísmica de Ambato viene dada por sismos corticales; este comportamiento se repite 
para otros períodos y para PR de 975 y 2475 años. Por esto se decidió trabajar con 
sismos asociados a fallas geológicas locales, ya que estos van a causar más daño 
en las estructuras, por estar más próximos a las fallas. 


En la actualidad, los análisis de peligro sísmico, tanto probabilístico como 
deterministas, requieren una evaluación del potencial de movimiento de origen 
tectónico futuro en la región ya que es cada vez más frecuente la inclusión de un 
modelo de fallas activas. Seleccionando todas las fuentes que pueden generar 
futuros sismos. 


El hecho de conocer directamente las fallas como fuentes sísmicas, permite 
incorporar nuevos parámetros al cálculo probabilista de la peligrosidad, como los 
mecanismos focales, geometría del plano de ruptura, tasas de deslizamiento, 
periodos de recurrencia, etc. Se analizará más adelante 


De acuerdo con los principios geológicos que se están desarrollado, el 
planteamiento que se sigue al modelizar las fallas como fuentes sísmicas 
independientes, en lugar de utilizar zonas sismogénicas, incide en dos aspectos: El 
cambio de la geometría de la fuente sísmica es decir las fallas entran en el cálculo 
como elementos independientes, ya sea como áreas que define el plano de falla o 
como líneas que definen la traza de la falla proyectada en la superficie y en segundo 
lugar el cambio de la fuente de información para caracterizar el potencial sísmico , 
la caracterización del potencial sísmico se efectúa a través de la tasa de ocurrencia 
de sismos N(m) y la proporción de estos en funcion de su tamaño £, por lo que la 
fuente básica es el catalogo sísmico, en el caso de las fallas los parámetros que 
caracterizan el potencial sísmico son difícilmente extraíbles del catálogo sísmico 
para cada falla individual y se derivan de las características físicas de la falla. (Rivas, 
2014) 


Al considerarse una falla como fuente sísmica activa e independiente, se 
modeliza el plano de la falla en el que se asume que se produce la ruptura, sabiendo 
que se generaran superficies irregulares de ruptura, al no poseer de forma precisa 
los datos geométricos de la falla de manera general se utiliza un plano rectangular. 
Los parámetros geométricos que se utilizan para tal modelización se muestran en 
la figura 4. 
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Figura 4 Geometria y parametros de la falla. Fuente: Rivas (2014). 


Donde: L es el largo de la falla; W es el ancho; D es el desplazamiento co- 
sísmico; 6 el azimut que varía entre O y 360°; 6 el angulo de buzamiento que varia 
entre O y 90%, también conocido como dip; A es el vector de deslizamiento que se 
encuentra entre 0 y 90° 


La ciudad de Ambato ubicada en la provincia de Tungurahua, en una zona 
volcanica y sismicamente activa, para esto vasta ver la figura 5, en que se presentan 
las fallas geologicas que atraviesan o estan muy cerca de area urbana y rural de la 
ciudad de Ambato, con circulo negro se indican las fallas capaces que pueden 
generar sismos de magnitud > 6 y con círculo rojo las que pueden producir un sismo 
de magnitud > 6.5. De estas, tres son las más preocupantes (para evaluar la 
peligrosidad sísmica en el área urbana) y son las mostradas en la figura 6; son las 
fallas de: Ambato, Huachi y Totora. Ahora, para evaluar la peligrosidad sísmica en 
el área urbana y rural de la ciudad de Ambato, se debe considerar las 7 fallas que 
están indicadas en la figura 7 





Figura 5 Fallas geológicas en radios de 20 y 50 km, que pueden afectar al área 
urbana y rural de la ciudad de Ambato. Los nudos capaces de generar 
terremotos con M26 son indicados con el círculo negro, mientras que los 
nudos capaces de generar terremotos con M26.5 son indicados con el 
círculo rojo. Fuente: Chunga y Mulas (2018) 
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El área rural de la ciudad de Ambato es por lo menos unas cuatro veces 
que el área urbana, por este motivo no se puede evaluar la peligrosidad sísmica 
solo con las fallas mostradas en la figura 6. Esas están muy bien para solo el área 
urbana (Aguiar y Rivas, 2018). 






















ka T pry A PERS 
y dS OIR RTS A j a ON 
A f i A J Y - 
de i A ac 
Whe nn EY 
ps f + + 
eo. AA i 
y qe di Falla tectónica: EC-344a F 
dé | Y) NOMBRE Huachi 
ir, Falla tectónica: EC-344b : 
i NG A tipo inversa ciega E 
$ NOMBRE Ambito i 
€ a i 
‘ 15, y Spa inversa ciega y Sa nta Rosa Longitud [km] 12,04 1 A 
4 13 yy E | Longitud [km] 16.30 af 7 \/ -P 
q Se ry Zoom to eee de 
40 a AL ; 
i að? > 4 y ) ` a 
Y? reía Mi 
~ er f 
Y A F 
f 167 
P oY + j 
) : ie 
Ý : 
~L £30 | 





Y No 
rt Quebrady Worejon 1 Nee 
i . 
e 


i Bartolome E 
le Pan > l 


Falla tectónica: Ec-345 


f 

| 
NOMBRE Totoras E 
tipo transcurrente inversa ciega 
Longitud [km] 17.65 | 


| Zoomto ... 

K \ 4 SA N 

Figura 6 Fallas geológicas seleccionadas en el estudio de peligrosidad sísmica del 
área urbana de Ambato. Fuente: Chunga y Mulas (2018); Eguez (2003). 


9860000 








30000 9840000 


476 Roberto Aguiar y Paola Serrano 


Figura 7 Fallas geologicas seleccionadas en el estudio para evaluar la peligrosidad 
sismica del area urbana y rural de Ambato. 


A mas de los resultados encontrados en la desagregacion de la peligrosidad 
sismica (figura 3) vale la pena recordar que el terremoto de Pedernales del 16 de 
abril de 2016 (M=7.8) en la ciudad de Ambato fue sentido como un sismo de baja 
intensidad que no causó daño en las estructuras. En cambio, si se registra un sismo 
de magnitud 6, en una de las fallas ciegas de la ciudad, el daño será cuantioso. 


En la tabla 1, se presentan las fallas que van a considerarse en el Método 
Determinístico para hallar los espectros ante la componente horizontal de 
movimiento del suelo. Ahí se indica la magnitud máxima esperada en cada falla 
hallada con las ecuaciones (1) y (2) propuestas por Leonard (2010) y la ecuación 
(3) de Wells & Coppersmith (1994); para calcular la magnitud se ha utilizado como 
parámetro el área (4) y la longitud de rotura en la superficie de la falla (L). Se 
encontró la magnitud con cada una de las ecuaciones, se obtuvo un promedio y se 
redondeó al valor inmediatamente superior. 


My = 1.52 log(L) + 4.4 (1) 
My = 1.0log(A) + 4.0 (2) 
M,, = 1.22 log(L) + 5.0 (3) 


Tabla 1 Parámetros de las Fallas ciegas que atraviesan el área urbana y rural de 
la ciudad de Ambato. 


A land 
(Km) | (Km?) | (Km) 


Ambato | 45 | 90 | 6.5 | 16.3 | 97.80 | 6 | Inversa Ciega | 
| Huachi | 45 | 90 | 6.3 | 12.04 | 72.24 | 6 | Inversa Ciega — 


Totoras 45 6.5 | 17.65 | 105.90 EE $ 
Inversa Ciega 


— F40 | 45 | 90 | 61/150 | 120 | 8 |Inversa_ | 

| F11 | 45 | 90 |66/ 110 | 99 | 9 [Inversa 
— F16 | 85 | 5 | 62] 180 | 144 | 8 | Cizalla Dextral | 
| Fee | 45 | 90160] 130 | 91 | 7 [Norma | 





3. CALCULO DEL PERIODO DE RECURRENCIA 


Ahora bien, interesa conocer el periodo de recurrencia de estos sismos de 
magnitud máxima y para ello se tienen dos formas, a saber: i) Método del Terremoto 
Característico; ii) Método de Gutenberg y Richter modificado, los mismos que se 
exponen a continuación. 


3.1 Método del Terremoto Característico 


El modelo de terremoto característico (TC) supone que en la falla solamente 
se producen rupturas completas de todo el plano de falla o segmento, por tanto, la 
falla solo rompe en forma de sismos de la misma magnitud, siendo esta la magnitud 
máxima. El modelo TC es un modelo simple que iguala el momento sísmico 
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acumulado en la falla durante un tiempo de recurrencia al momento sísmico liberado 
en el máximo terremoto de la falla. (Rivas, 2014) 

Por medio de la expresión de Hanks y Kanamori (1979), se puede 
establecer una relación entre el momento sísmico y la magnitud máxima de dicho 
terremoto expresada en magnitud momento Mw, la cual permite conocer la magnitud 
máxima que se espera en la fuente y la energía que se libera en el terremoto. 


(My +10.7)3 


: a (4) 


Haciendo una suposicion de que en todo el plano de falla se esta 
acumulando energia uniformemente, la tasa de momento sismico M, que se 
acumula en la falla está relacionada con la tasa de deslizamiento anual de la falla ù 
según la ecuación de Brune (1968). 


kgf 


a u es la tasa 


Donde u es el módulo de rigidez o de cizalla (u ~ 3 * 10? 
de deslizamiento y 4 es el área del plano de falla. 


El modelo TC es un modelo totalmente periódico, en el que la energía 
acumulada es igual a la energía libertad en el TC. Por lo tanto, conociendo el Mọ 
que se liberaría en el terremoto y la tasa de acumulación anual de momento sísmico 
de la falla Mọ, se puede conocer el período de recurrencia Tr, es decir el tiempo 
que transcurre entre dos eventos característicos con la siguiente ecuación. 


Mo 
Mo 

En la tabla 2 se presenta el Periodo de Recurrencia hallado con el modelo 

del Terremoto Característico, para las fallas que se han indicado en la tabla 1 y 


figura 7, trabajando con una tasa de movimiento de 3.5 mm/año y un módulo de 


| k 
deformación en roca de 3 * 107? “or 


Tabla 2 Períodos de recurrencia encontrados con el Método del Terremoto 
Característico para las fallas consideradas en el estudio. 


Mo Mo A=L*w Mo Mo T,. 
(Dyn*cm) | (N*m) (m?) (kgf m/afio) | (Nm/año) | (años) 
3.16*1025 | 3.16*10%8 | 72.24* 10% | 7.585*1015 | 7.585*10* 


6.3*1025 | 6.3*1018 | 97.8*10° | 1.027*1015 | 1.027*101° | 626.5 
6.3*1025 | 6.3*1018 | 105.9*10° | 1.112*1075 | 1.112*10° | 578.6 


1.512% 101° 
F28 





3.2 Modelo de Gutenberg y Richter modificado 
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Con el propósito de determinar períodos de recurrencia, para una falla 
geológica, otros autores contemplan a la utilización de otros modelos como es el de 
Gutenberg y Richter modificado, la cual establece una proporción constante entre 
el número de sismos grandes y pequeños para una zona y en un periodo, con la 
posibilidad de que la falla produzca terremotos de diferente magnitud con 
periodos de ocurrencia diferentes. 


Siendo formulada por Charles Francis Ritchter y Beno Gutenberg en 1944 
y modificadas por Consentino et al. (1977), en la cual se integró parámetros como; 
la tasa acumulada de sismos (N(,,,), la tasa acumulada de sismos de magnitud 


mayor o igual a Min (Nein) mostrada en la siguiente ecuacion: 


, , e B(m) = eB (max) 
Nom) > NM min g e ~B(M min) — e"B(Mmax) (7) 


4 < __Mal@= A) (oP Hn — 8 Mmas) 
vee f [e-f Mmax Mo (Mmax) — e=B Mmin Mo(Mmin)! 


Mo(Mmin) = e(Ctd Mmin) M)(Mmax) = e(c+d Mmax) (9) 


č = 16.05 x In(10) d = 1.5 » ln (10) (10) 


La tasa acumulada de momento sísmico será liberada por medio de sismos 
de diferente magnitud, las cuales varían desde una magnitud mínima Mmin, hasta 
una magnitud máxima M max, la cual viene definida por el tamaño del plano de falla. 
Aguiar et al. (2015). 


Las fallas indicadas en la figura 7, se encuentran en la zona fuente 6, 
conocida como Falla de Pallatanga, del modelo de Beauval et al. (2018). Se hace 
referencia a este trabajo ya que se va a trabajar con los parámetros de sismicidad 
a = 2.8012; b = 0.73, del modelo clásico de Gutenberg y Richter log N = a — b M. 
Los valores anotados fueron hallados para la zona fuente 6. 


Tabla 3 Período de recurrencia, hallada con el modelo de GR modificado en la falla 
Ambato. 


TASA ANUAL NO 
ACUMULADA PERÍODO DE 


RECURRENCIA 
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En la tabla 3, se presenta el calculo del periodo de recurrencia, en la falla 
Ambato y a continuacion se describe su contenido, para entender la forma de 
calculo. En la primera columna se presenta la magnitud (m) considerada; en la 
segunda columna la tasa anual acumulada N(m) de sismos mayor a la magnitud 
(m); en la tercera columna, la tasa anual Nan) no acumulada de sismos que se 
obtiene restando la tasa anual acumulada de dos magnitudes consecutivas, 
finalmente se tiene el período de recurrencia que no es más que la inversa de la 
tasa anual no acumulada. 


Tabla 4 Períodos de recurrencia, hallado con modelo de GR modificado para las 
fallas de Huachi y Totora. 


Falla de Huachi Falla de Totoras 


E a e a 


| 6.0 | 0.0017 | 0.0017 | 276.43 | 6.0 | 0.0025 | 0.0025 | 398.73 | 





Tabla 5 Períodos de recurrencia, hallado con modelo de GR modificado para las 
fallas F10 y F11. 


Falla F100 Fla Fd S 
om | Nm | Nom | r | m | Nm | Nom | T 
60 | 0.0030 | 0.0030 | 338.63 | 6.0 | 0.0023 | 0.0023 | 201.40 | 
64 | 0.0000 | 0.0000 | 551.44 | 6.3 | 0.0000 | 0.0000 | 274.40 | 





Tabla 6 Períodos de recurrencia, hallado con modelo de GR modificado para las 
fallas F16 y F28. 


o Falae | Fallar2s ë O—ë ë | 
Oom | Nm | Nom | r | m | Nom | Nom | r 


60 | 0.0034 | 0.0034 | 293.51 | 6.0 | 0.0022 | 0.0022 | 216.89 | 





En las tablas 4 a 6 se presenta el período de retorno para las otras 6 fallas. 
El período de recurrencia para tener un sismo con magnitud máxima en alguna de 


480 Roberto Aguiar y Paola Serrano 


las fallas que pueden afectar al area urbana y rural de Ambato varia entre 274 y 676 
anos. 


4. ECUACIONES DE MOVIMIENTOS FUERTES 


Las ecuaciones de predicción de movimientos fuertes, en español o Ground 
Motion Prediction Equations GMPEs, son fundamentales para la realización de 
estudios de peligrosidad sísmica, es así que día a día aparecen nuevas 
contribuciones científicas para tratar de predecir la atenuación del movimiento del 
suelo con la distancia, es complejo el tema ya que hay demasiadas variables 
involucradas en el cálculo. 


Estas ecuaciones son semi empíricas que tienen un gran conocimiento de 
propagación de ondas y una gran cantidad de registros sísmicos, por la gran 
cantidad de terremotos que se registran en el mundo, un nuevo sismo aporta con 
una gran cantidad de registros de aceleraciones, velocidades y desplazamientos del 
suelo, en las diferentes estaciones donde se encuentran los equipos. Información 
que para los científicos es fundamental para mejorar sus modelos de predicción. En 
Douglas (2019) se presenta un resumen muy corto de cada modelo publicado entre 
1964 y 2019, es un libro de 641 páginas. 


De tal manera, que la primera actividad que se debe hacer es imponerse 
criterios, para seleccionar de entre tantos modelos de movimientos fuertes, cuáles 
son los más adecuados teniendo en cuenta la realidad sísmica del Ecuador, el tipo 
de suelo que se tiene en la sierra, por ejemplo, predomina la “cangahua” entonces 
se deberán buscar modelos cuya base de datos sean en ese tipo de suelo. 


Para el estudio de peligrosidad sísmica del área urbana de la ciudad de 
Ambato, se hizo este estudio (Aguiar y Rivas, 2018) y se seleccionaron los 
siguientes modelos: Campbell y Bozorgnia (2014), que se abreviara C&B; 
Abrahamson, Silva y Kamai (2014) ASK; Choiu y Youngs (2014) CY. 


Ahora, a esta lista para el estudio de la peligrosidad sísmica del área urbana 
y rural de la ciudad de Ambato se añaden los modelos de Zhao et al. (2016) porque 
buena parte de la base de datos son en cangahua. Este es el modelo más difícil de 
los 6 que se van a utilizar y es un contraste con el modelo muy elemental de Zhao 
et al. (2006). En el nuevo modelo incorporan entre sus variables el vulcanismo, esto 
lo hace más interesante ya que es uno de los pocos modelos que existen en el 
mundo que consideran a los volcanes. 


La lista de modelos de movimientos fuertes se complementa con los 
trabajos de Boore y Atkinson (2008) B8A, y el modelo de Idriss (2014) que se 
llamará con el nombre del autor. 


Por lo tanto, se consideran 6 modelos para encontrar las aceleraciones 
espectrales en el área urbana y rural de Ambato, empleando la metodología 
determinística. 


Los tres primeros modelos están detallados en Aguiar y Rivas (2018); 
Aguiar (2016) por lo que se decidió en este artículo presentar sus ecuaciones de 
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calculo, sin entrar en mayor detalle y presentar tablas-comparativas, que ayuden a 
tener una idea de estos modelos. 


e Modelo de: Campbell y Bozorgnia (2014) 


InY = Imag + fais + frit T Fang + fsite + fsea + Tiso + fain + fatn (11) 


Donde Y es la aceleración espectral expresada como una fracción de la 
gravedad; fray es una función de la magnitud M,,; fa;, es una función de la magnitud 
y de la distancia del sitio al plano de ruptura R,,,,,; ff, depende de la magnitud y del 
tipo de falla; fang es una función que depende del efecto hanginwall; f.¡¿¿ es una 
función que toma en cuenta los efectos de sitio y depende fundamentalmente de 
V¿3y velocidad de la onda de corte en los primeros 30 metros; f... es una función 
que toma en cuenta el efecto de cuenca, depende de Z, que es la profundidad 
donde la velocidad de la onda de corte es 2.5 km/S; fhyp es una función que depende 
de la profundidad focal Zhyp; faip es una función que depende del ángulo de 
buzamiento de la falla 6 y de la magnitud; f.,.,, es una función de la atenuación 
anelástica. 


e Modelo de: Abrahamson, Silva y Kamai (2014) 
ln Y = 
fı (M, Rrup) + Frv: f7 (M) + Fn: fe (M) + Fas: f11(CRjp) 
+ fs (Sa1100, Vs30) + Fuw: f4 (Rjb, Rrup, Rx, Ryo, W, dip, Zror, M) + fe (Zror) 
+ fio (Z1.0, Vs30) + Regional (Vs30, Rrup) (12) 


Donde f} es una función de magnitud y distancia; f, y f¿ son funciones del 
tipo de falla; f,, es una función de réplicas; f es una función de sitio en función de 
Vozo; f4 es una función de hanginwall; f£¿ es una función de Zror que es la 
profundidad al borde superior del plano de ruptura; f, es una función del efecto de 
cuenca; Regional es una función que permite trabajar mediante variables 
indicadoras con datos de determinadas regiones del mundo de donde proviene la 
información sísmica. 


e Modelo de Choiu y Youngs (2014) 


La mayor parte de funciones que intervienen en su evaluacion son similares 
a la de los modelos anteriores. Un aspecto fundamental de este modelo y que no 
consideran los modelos anteriores es la directividad de las ondas sismicas en la 
función Appp. 


Este modelo a más de Zo», considera la función AZ ror que ayuda a definir 
de mejor forma al primer valor considerando el valor medio de esta profundidad en 
función del tipo de falla y de la magnitud. 


A continuación, se indica la ecuación de cálculo para hallar Y, que es la 
aceleración en roca expresada como una fracción de la gravedad. 
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In(Yreri) = C1 + (Ca + el Fryi 
+ fcin + aoe ee Fumi 
F fe, + ea a AZ Tori 
+ z ag E a. (cos 6;)* 


C2 — C3 





+C2(M; — 6) + In(1 + eon'em-Mi ) 


n 


+ Ca n(Rrori + cs cosh(c¿ * max(M; — Cym, 0))) 


Cy2 
+ (C4a — C In ( R} + GR) + fe + j 
(csa = c4) de RE 2" cosh(max(M — Cy3,0)) see 


+ Cg max|{ 1 


max(Reypi; — 40,0) =~ fmax(M; — 5.5, 0) E _ 
Rea Rīsi t ZīoRi 
ae 5; f 1— coa )t (=) ie ooo 
+ Co FywijCOS Ôi }Coa + (1 — Coa tan a Read 
(13) 
e Modelo de Boore y Atkinson (2008) 
In Y = Fy(M) + Fp(Ryg, M) + Es(Vozo,Ry8, M) (14) 


Donde Fy es una función de magnitud; F, es la función de distancia; F, es 
la función de sitio. En este modelo al igual que en los anteriores se trabaja con R;g 
es la distancia Joyner-Boore, definida como la distancia más cercana a la 
proyección en la superficie de la falla, la cual es aproximadamente igual a la 
distancia epicentral para eventos de magnitud seis. 


e Modelo de Idriss (2014) 
In(Y) = a, + a,(M) + az(8.5 — My? — (8, + B2M) * Lo(Ryop + 10) + eLn(Vgz0 ) + Y Rrup + GF 


(15) 


En forma explícita no se indican las funciones, pero se ve que es función de 
la magnitud; de la distancia (Ras). de los efectos de sitio con el parámetro (V.z. ) y 
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el tipo de falla mediante la variable indicadora F que puede valer: O (Falla 
transcurrente); 1 (Falla Normal); 2 (Falla Inversa); 3 (Inversa oblicua); 4 (Normal 
oblicua); en este sentido este modelo permite considerar 5 tipos de fallas, es más 
específico que los otros en cuanto se refiere al tipo de falla. 


e Modelo de Zhao et al. (2016) 


El modelo presenta dos ecuaciones una para sismos de corteza (fuente 
menor a 25 km) y otra para sismos de manto superior (fuente mayor a 25 km) 


Sismos de corteza 


log. (y; ;) = fier + Ger log. (ri ;) 
+ Gert log. (xi; + 200.0) + gn(xi;) Ped Cor Xi y + Ver 
+ log.(A™) + ij + ni (16) 


Sismos de manto superior 


loge(yi,;) = frum + Jum loge (rij) 
+ geri log. (xi; + 200.0) + gn(xi;) eye Cor Xi y + Yer 
Flog AA IFs EN (17) 


El subíndice cr indica que se trata de sismos de corteza y el subíndice um 
de manto superior, con esta aclaración se indican las funciones sin el subíndice; fn 
representa el efecto de fuente que es función de la magnitud, de la profundidad focal 
y del tipo de falla; ger, gum es la tasa de atenuación geométrica que es función de 
la distancia y la magnitud, el modelo considera ecuaciones para el campo cercano 
(gn) y para larga distancia (ger; ); ecr, €um, ES la tasa de atenuación anelastica; e% 
es la tasa de atenuación anelástica asociada con la ruta volcánica (esto se comenta 


mas adelante); A®’, A“, son funciones de sitio. 


Por ser este modelo más reciente, en comparación con los otros que se han 
presentado en este apartado, vale la pena mencionar dos aspectos que se 
consideran importantes, ellos son: 


La clasificación de los suelos, que ellos consideran en su estudio, que se 
indica en la tabla 7. Al margen del tema que se está tratando, es importante ver que, 
en la clasificación de los perfiles de suelo sísmico a más de la velocidad de la onda 
de corte, incorporan el período del suelo como variable fundamental. Esto es 
fundamental ya que se puede tener un suelo que de acuerdo al V,3, clasifique como 
tipo C (NEC-15) que es bueno, pero si tiene periodos altos, va a ver amplificacion 
de las ondas sismicas. 


Tabla 7 Clasificación de suelos considerada por Zhao et al. (2015) y similitud con 
NEHRP 
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Perfil de [Descripción — | =. de suelo 
suelo (m/s) NEHRP 


C 


SC Il Suelo Duro 0.2<T<04s 300 < Viz 
< 600 


SC Ill Suelo 0.4 <T < 0.65s 200 < Vozo 
Mediano < 300 


SC IV Va < 200 


SC IV1 0.6<T<10s 120 < Vozo E 
< 200 


SCIv2 | | T>108 |  Vsp<120 





En la puntilla de Bahía de Caráquez en Ecuador se dio que los suelos, por 
el valor del V,. son tipo C (clasificación de NEC-15), pero tienen períodos de suelo 
altos alrededor de 0.9 s, lo que ocasionó amplificación de las ondas sísmicas y con 
ello se tuvo gran daño en las estructuras durante el terremoto de Pedernales de 
2016. Aguiar (2017) 


Otro aspecto importante del modelo de Zhao et al. (2016) tiene que ver con 
la consideración del vulcanismo en el modelo, para ello se definen 4 casos, que se 
muestran en la figura 8, y en forma resumida son: i) Epicentro fuera de la zona 
volcánica (ZV) y estación o punto de interés dentro de ZV; ii) Estación y epicentro 
dentro de ZV; iii) Epicentro en ZV y estación fuera de ZV; iv) Estación y epicentro 
fuera de ZV pero que atraviesa la ZV. En todos los casos la distancia x” queda 
dentro de ZV. 







Station 


Ground 


Earthquake 


Figura 8 a) Definición de la trayectoria volcánica para cuatro casos; b) Distancia 
volcánica horizontal e inclinada. 


En la figura 7, se ve que las parroquias urbanas y rurales de la 
ciudad de Ambato se hallan bastante distantes del Volcán Tungurahua, con 
lo que la Zona Volcánica no llega a la ciudad. En caso de una erupción lo 
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que llega es ceniza volcanica. Sin embargo, de esto el modelo recomienda que 
como mínimo x” = 12 km., cantidad que se consideró en el cálculo para hallar 
espectros. 


Ahora, con toda la información indicada en el presente apartado, en la tabla 
8 se indican las variables involucradas en su formulación con el propósito de ver 
cuál es la más completa. Para no alargar el texto, en la columna uno, se tiene el 
parámetro M se debe interpretar como que, si el modelo considera la función de 
magnitud, en este caso todos lo consideran. 


Tabla 8 Parámetros básicos que consideran los modelos seleccionados para el 
estudio. 


C&B-14 | C&Y-14 AAA a A: 
a 


| Zrop | X | 
| Fry Fy, SS |X 


Ryu O 


X| X| X< 


an 


O X 
Oa 


O Rap | X 
O R | X | 
o Roe | = 


Zas OA] XX | XX | XX | ~ | XX 
Xx 


Zhyp 





El significado de los parámetros indicados en la tabla 8, es el siguiente: M 
es la función de magnitud; Z y es la distancia desde la superficie al borde superior 
de la falla; F,y función de falla inversa; Fyy es la función de falla normal; SS es la 
función de falla transcurrente; R,.,,, es la distancia del sitio de interés al plano de 
falla o la distancia más corta a la falla; R; es la distancia en el plano horizontal 
desde el sitio a la parte inferior de la falla, ver figura 8; R, es la distancia horizontal 
desde el sitio a la parte superior de la falla; R,,, es la distancia al final de la ruptura, 
paralela a la falla; F,,y es la función de Hanging Wall; V.-. es la velocidad de la onda 
de corte en los 30 primeros metros; Z, s, Z,, es la profundidad a la cual la velocidad 
de la onda de corte es 2500 m/s, o 1000 m/s; Z;,,y es la profundidad al hipocentro; 
X” es la distancia entre el epicentro y el sitio en la Zona Volcánica. 


Se recomienda mirar con detenimiento la figura 9 en que se presentan 
algunas de las variables involucradas en las ecuaciones de movimientos fuertes que 
se han presentado. 


Por otra parte, en el Anexo 1, se indica la base de datos, consideradas en 
cada modelo. 


486 Roberto Aguiar y Paola Serrano 





'——— (0) A p0) aaeaaaeaeaeaeututliet 
Sto i Rd — m 





Figura 9 Algunas de las variables involucradas en los modelos de movimientos 
fuertes. 


5. MALLA DE PUNTOS CONSIDERADA EN EL ESTUDIO 


En la microzonificación sísmica del área urbana de la ciudad de Ambato se 
clasificaron los suelos de acuerdo al NEC-15, por ello los puntos que se presentan 
en la figura 10, con color verde-amarillo corresponden al perfil de suelos C. (Aguiar 
y Rivas, 2018) 


MALLA DE PUNTOS AMBATO 
9885000 
9880000 
9875000 
9870000 
9865000 
9860000 
9855000 
9850000 
9845000 


5840000 





9835000 


9830000 


720000 730000 740000 750000 760000 770000 780000 


© ZONA RURAL ZONA URBANA 


Figura 10 Malla considerada en el estudio; con color azul se indica el area rural y 
con verde amarillo el area urbana de la ciudad de Ambato. 


Para el area rural, no se han realizado estudios geofisicos y geotecnicos 
que permitan clasificar los suelos, por lo que se asume que todos son perfiles de 
suelo C. Con esta acotacion, en la figura 10 se presenta la malla de puntos en las 
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que se obtuvieron espectros de respuesta para 5% de amortiguamiento utilizando 
las ecuaciones de movimientos fuertes indicadas en el apartado anterior. 


La separacion de los puntos de la malla es de 500 m, en los dos sentidos, 
se encontró espectros en cada punto, considerando que un sismo de magnitud 
máxima (ver tabla 1) luego se obtuvo un promedio para tres casos: i) Considerando 
solo los puntos del área urbana; li) Lo propio, pero del área rural; y iii) Todos los 
puntos de la malla; espectros que se presentan en el siguiente apartado. 


6. RESULTADOS 


En la figura 11 a, se presentan los espectros promedios, encontrados con 
cada uno de los modelos seleccionados en el área urbana, ante un sismo de 
magnitud 6.5, en la falla Ambato y en la figura 11 b, se tiene los espectros promedios 
para el área rural. Al estar la falla Ambato más cerca del área urbana las ordenadas 
espectrales son mayores. 
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Figura 11 Espectros para un sismo de magnitud 6.5 en la falla Ambato; a) Area 
Urbana; b) Area Rural. 
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En la figura 11, se ha presentado espectros medios, cuyo nivel de confianza 
es del 54% y espectros medios mas una desviacion estandar que tiene una 
confiabilidad del 84%. Por cierto, en las ecuaciones 11 a 17, no se colocó la 
incertidumbre de los modelos que se los cubre con la desviación estándar. 


En caso de generarse un sismo de magnitud 6.3, en la falla Huachi, los 
espectros promedios que se esperan en el área urbana y rural de la ciudad se 
indican en la figura 12 a, y b, respectivamente. 


Nuevamente por la cercanía de la falla Huachi al área urbana los espectros 
son mayores (figura 12 a) que los del área rural (figura 12 b) Esto significa que habrá 
mayor daño en el área urbana que en la rural, ante un sismo asociado a la falla 
Huachi. 
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Figura 12 Espectros para un sismo de magnitud 6.3 en la falla Huachi; a) Area 
Urbana; b) Area Rural. 


Ahora, en la figura 13, se presentan los espectros de respuesta elásticos 
para el caso de que se produzca un sismo de magnitud 6.5; los espectros de la 
figura 13 a, son para el área urbana y los de la 13 b, para el área rural. 
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Figura 13 Espectros para un sismo de magnitud 6.5 en la falla Totora; a) Area 
Urbana; b) Area Rural. 


Las tres fallas: Ambato, Huachi y Totora, se consideraron en el estudio de 
microzonificación sísmica del área urbana en 2018, y se ve en las figuras 11 a 13, 
que estas fallas no generan espectros altos en el área rural. 


Ahora se van a considerar otras fallas, algunas de ellas se encuentran en el 
limite del área urbana como son las fallas F10 y F11. En la figura 14 se presentan 
los espectros asociados a la Falla F10; como se encuentra en el límite entre el área 
urbana y rural sus espectros son parecidos en estas dos áreas y tiene menores 
valores a las encontradas con las tres fallas precedentes, de tal manera que estuvo 
bien no considerarla en el estudio realizado en 2018, para el área urbana. 
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Figura 14 Espectros para un sismo de magnitud 6.1 en la falla F10; a) Area Urbana; 
b) Area Rural. 


La falla F11, si presenta valores altos de los espectros que son comparables 
con los espectros de las fallas de Ambato, Huachi y Totora. Esto no se esperaba y 
vale la pena actualizar el estudio de microzonificación sísmica del área urbana de 
Ambato, habida cuenta que se dispone de una clasificación de los suelos. 


Para el área rural, no se dispone de estudios geofísicos y geotécnicos que 
permitan saber con detalle los sitios donde se tienen suelos C, D, E. En este artículo 
se ha considerado que toda el área rural es tipo C, porque el objetivo es ver cuál es 
el comportamiento de otros modelos de movimientos fuertes y la gran expectativa 
que se tenía es el modelo de Zhao et al. (2016), que dicho sea de paso es el modelo 
más complejo de los 6 utilizados en el estudio. 


Para valores medios, confiabilidad del 50% (sin desviación estándar) la 
ecuación de Zhao et al. (2016) reporta ordenadas espectrales que se encuentran 
por la media a las que se hallan con los otros modelos, tal vez un poco menos de la 
media, pero para confiabilidad de 84% (con desviación estándar) están un poco más 
arriba de la media. 
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Figura 15 Espectros para un sismo de magnitud 6.6 en la falla F11; a) Área Urbana; 
b) Area Rural. 


Los espectros que se generan con las diferentes ecuaciones de 
movimientos fuertes ante un sismo en la falla F16, de magnitud 6.2, se observan en 
la figura 16 a, para el área urbana y en la figura 16 b, para el área rural. Al comparar 


estos grupos de espectros se ve que los del área rural, son ligeramente mayores a 
los del área urbana. 


Los espectros del área urbana hallados con un sismo en la falla F16 son 
menores a los obtenidos en la microzonificación sísmica de Ambato de 2018. 


Tanto para el área urbana como rural, los espectros más altos se hallan con 
el modelo de Idris (2014) este comportamiento también se lo ha observado en los 
espectros generados por sismos en las otras fallas consideradas en el estudio. 


Para el área rural los espectros obtenidos ante un sismo de magnitud 6.2 
en la falla F16 son los más altos, con las 7 fallas consideradas en el estudio. 
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Figura 16 Espectros para un sismo de magnitud 6.2 en la falla F16; a) Area Urbana; 
b) Area Rural. 


Resta por indicar los espectros promedios que se obtienen a partir de un 
sismo de magnitud 6, en la falla F28. En la figura 17 a, se ven los espectros para el 
area urbana, los mismos que son comparables a los encontrados en la 
microzonificación sísmica de 2018, para el area urbana. 


Para el área rural, son bajos los espectros generados por un hipotético 
sismo en la falla F28. 


En todas las fallas se consideró que el hipocentro de los sismos se da en el 
punto medio del volumen de litósfera de la falla. Se encontró los espectros en cada 
punto de las mallas, para la magnitud máxima esperada en cada falla y después se 
obtuvieron los espectros promedios, de tal manera que es un trabajo largo. 


Se utilizó un programa desarrollado en EXCEL por varios autores, la primera 
de ellas fue Rivas (2014); después Trujillo (2015) lo adaptó a las fallas ciegas de 
Quito; luego Quizanga (2015) como parte de su Maestría en Ciencias, hizo 
significativos aportes todo a las fallas de Quito; nuevamente Quizanga (2018) 
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trabajo con este programa, pero esta vez adaptado a las fallas de Ambato, Huachi 
y Totoras. Finalmente, Serrano (2020) ha incorporado las fallas F10, F11, F16 y 
F28; al igual que las ecuaciones de movimientos fuertes de Boore y Atkinson (2008), 
Zhao et al. (2016) e Idriss (2014). 


En todo este largo recorrido de la evolución del programa, acorde al 
desarrollo científico, el primero de los autores de este artículo ha estado 
supervisando el desarrollo y analizando los nuevos modelos a utilizar. De tal manera 
que el programa tiene un sustento bastante sólido. 
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Figura 17 Espectros para un sismo de magnitud 6.0 en la falla F28; a) Área Urbana; 
b) Area Rural. 
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7. COMENTARIOS Y CONCLUSIONES 


En el 2018, se realizó la microzonificación sísmica del área urbana de la 
ciudad de Ambato, para lo cual se conformó un equipo de trabajo multidisciplinario 
y se hizo un gran trabajo de campo en el área de la geofísica, geotécnica y la 
geología; este trabajo sirvió de base para la clasificación sísmica de los suelos y 
para los estudios de peligrosidad sísmica en forma probabilística y determinística 
para encontrar espectros de diseño y máximo considerado de: aceleraciones para 
la componente horizontal y vertical del movimiento del suelo, así como para obtener 
espectros de desplazamiento. Aguiar y Rivas (2018) 


Ahora, se retomó el tema en forma exploratoria por decirlo de alguna 
manera incorporando el área rural de la ciudad de Ambato, que prácticamente es 
cuatro veces el área urbana y se encuentra muy poblada. 


Es exploratoria porque no se han realizado estudios preliminares básicos 
indicados en el primer párrafo de este apartado, por lo que se decidió suponer que 
toda el área rural es un perfil de suelo C; lo que interesaba es mostrar la metodología 
de cálculo de la peligrosidad sísmica en forma determinística, incorporar nuevas 
fallas al modelo, determinar los períodos de recurrencia con dos métodos a saber 
el del terremoto característico y el de Gutenberg y Richter modificado. 


Se ha incorporado nuevas ecuaciones de movimientos fuertes de tal 
manera que este momento se disponen de 6 ecuaciones, a pesar de que se 
obtuvieron resultados razonables se va a eliminar uno de los modelos e incorporar 
otro, para de esta forma ampliar en un futuro cercano la microzonificación sísmica 
del área urbana y rural de Ambato. 
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Anexo 1 Base de datos de los modelos considerados 


: Z et al. 
arámetro So r 4 SA 4 a aw | 


Magnitud Mínima lo 50 | 
Magnitud Máxima A a a 


Magnitud máxima 8.50 8.50 8.50 8.50 

falla transcurrente | 
een 8.00 8.00 8.50 8.50 8.00 | 
falla inversa | 
Magnitud máxima 7 50 8.00 7.00 8.00 

falla normal | 
Número de registros 8313 12 244 26 1658 70 | 
de fallas Normales | 
Número de registros 322 2.587 79 3117 1490 

de fallas Inversas | 


Número de registros 
de fallas 1561 221 1182 899 
transcurrentes 


California, Taiwan, 
California, California, J. USA. USA, Canada, 
Datos de Sismos wordwide | California | Japon, China, ia California Georgia, Greece, 
earthquakes Taiwan Iran, Italy Mexico, 
Turke 


(0-300) | (0-300) | (0-300) ES (0-400) | (02-175 
(150-1500) |(180-1500)| (180-1500) > 600 > 1300 Ta 
ZTOR 
eae Pett 
ee (afee | slee] 1 
(Km) 
(0-10) | (0-10) os [|x| 
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